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'^ '. Abstract 

cn ' In the standard model ( SM ) of electroweak interactions, CP noninvariance arises 
J^ ■ from the nonzero phase in the CKM matrix. Its contribution to the quark electric dipole 
t^^ ■ moment ( EDM ) vanishes surprisingly at two loop order. This makes the quark EDM 
^ ', extremely small in the SM. In this paper, we consider the two Higgs doublet extension 
(^ [ of the SM and assume that CP noninvariance is still encoded in the CKM matrix. We 
^ ] calculate the charged Higgs boson contribution to the quark EDM which naively should be 
■*-««. ■ of order eG|.(5(47r)~^m„(d)mfmfm^^ for the up ( down ) quark with possible enhancement 
p • factors of tan (3. Here 6 is the rephasing invariant of CP violation. However, contrary to 
A • the above naive expectation, we find that the charged Higgs boson contribution vanishes 
D ! strictly at two loop order. We show explicitly how this comes about and explains how it is 



^ 



related to the general form of Yukawa couplings in a spontaneously broken gauge theory. 
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I. INTRODUCTION 

CP noninvariance in any models with CPT symmetry will in general induce P- and 
T-violating electric dipole moments ( EDM ) for elementary particles through quantum 
effects. The discovery of these moments would be direct evidence of CP noninvariance 
outside the scope of the neutral kaon system and would help us identify the origin of CP 
noninvariance |l|. Of special interest among these moments is the neutron EDM. The 
current experimental upper bound for the neutron is |(i(A^)| < 1.1 x 10~^^ e cm and 
this limit is hopefully to be improved by several orders of magnitude in the near future . 

From theoretical point of view, any calculation of the neutron EDM proceeds in two 
steps. In the first step, one writes down all relevant operators which break the CP 
symmetry and involve only light degrees of freedom. These operators are usually defined 
at a high energy scale where CP violation occurs. They include the electric or chromo- 
electric dipole moments of the light quarks and gluons |^ , four quark operators ^ and 
possibly others. In the second step, they are evolved down to the typical hadronic scale 
and their effects on the neutron EDM are then calculated. While naive dimensional 
analysis or hadronic models have to be invoked in the second step, the first step can be 
implemented unambiguously once the model of CP violation is specified. In this work, 
we will be concerned with this first step calculation, in particular, the calculation of the 
quark EDM. 

In the standard model ( SM ) of electroweak interactions, CP noninvariance arises 
from the nonzero phase in the CKM matrix. The quark EDM vanishes trivially at one 
loop order since only moduli of the matrix are involved in the relevant amplitude. At 
two loop order, the fiavour structure is rich enough which could in principle allow for 
a CP- violating EDM 0; however, the final contribution to the quark EDM vanishes 
surprisingly when the sum over internal virtual fiavours is taken ||^ |^. This circumstance 
also appears in the W^ EDM 0] |l^. Although there are attempts to understand 
the vanishing result |]11| [|1^, it seems clear by now that it does not result from any 



symmetry which would dictate the zero EDM automatically at the lowest nontrivial order. 
Furthermore the vanishing result is accidental in some sense, as a result of specific Lorentz 
and fiavour structure. First, the W^ boson still acquires a nonvanishing, P- and T- 
violating magnetic quadrupole moment at two loop order |T^ though its EDM vanishes at 



the same order. Second, the quark EDM does not vanish any more when QCD corrections 
are included [Q [11| |15]. Of course, this makes the quark EDM extremely small in the 



SM. If the light quark EDM is one of the important contributions to the neutron EDM, it is 
then hopeless to observe the neutron EDM in the near future. Considering this, we would 
like to investigate how the quark EDM could be enhanced beyond the SM. In this paper 



we study the two Higgs doublet extension of the SM ^^. We will assume conservatively 
that CP noninvariance is still encoded in the CKM matrix so that the charged Higgs 
boson H^ is the only other particle besides W^ that mediates CP violation. Since the 
Yukawa couplings between H^ and quarks generally involve the relevant quark masses 
and are thus less universal as compared to the gauge couplings between W^ and quarks, 
we would expect naively that the H^ exchange will make a contribution to the quark 



EDM which should be of order d{u{d)) ~ eGp5{ATc) Tnu(d)'m't''^b''^H ^^ the up ( down 



quark. [ Here S is the rephasing invariant of CP violation |T^. ] As we will display later 
on, d{d) could even be enhanced by a large factor of tan^ (3 which would make d{d) easily 
reach the level of 10~^^ e cm for a charged Higgs mass of 200 GeV. However, a detailed 
calculation shows that the above naive expectation is actually not realized in the final 
result: the H^ contribution vanishes strictly at two loop order when the sum over internal 
virtual quark flavours is taken. We show explicitly how this null result comes about as 
a consequence of the general form of Yukawa couplings in a spontaneously broken gauge 
theory. 

The following sections are organized as follows. In section 2, we first discuss the renor- 
malization of one loop elements to be used in the complete two loop calculation. Then, we 
present respectively the contributions from exchanges of two charged Higgs bosons, one 
charged Higgs boson and one W^, for the general form of Yukawa couplings. In passing we 
also give the result from exchanges of two W^ which was previously calculated in the SM. 
The naive expectation for d{u) and d{d) is then verified. In section 3, we first show that 
the contribution involving exchanges of H^ or W^ vanishes when we sum over internal 
quark flavours. We also indicate a difference of the cancellation mechanism in the present 
case and in the SM case as computed in the unitarity gauge. Then we examine generally 
how it could be possible to have such a vanishing result. Section 4 is a recapitulation of 
our result. 



II. EXPLICIT RESULT OF CHARGED HIGGS BOSON CONTRIBUTIONS 

We shall calculate in this section the quark EDM arising from exchange of H^ and 
W^ in the two Higgs doublet extension of the SM. The Feynman diagrams at the lowest 
two loop order are depicted in Fig. 1. We shall denote the internal up- type quarks by 
Greek letters a, (3 etc, the internal down-type quarks by Latin letters j, k etc, and the 
external up- or down-type quark by e. Within this section, the flavours of internal quarks 
are fixed. We shall examine in the next section what will happen when summation over 
flavours is done. 

To set up our notation, we first list the relevant Feynman rules. The Feynman rule 
for the UadjW^ vertex is 

^^^ai7/.(l-75), (1) 

where g is the SU{2) weak coupling constant and Vaj is the entry (a,j) of the CKM 
matrix. Then the vertex djUaW~ is i — -/=V*j'y^{l — 75). The Feynman rule for the 



UadjH^ vertex is parametrized as 



^ V;,(C.,+C'-75), (2) 



2\/2mw 



where Caj and C'^^ are real constants and may depend on the masses of Ua and dj. The 

vertex djUaH~ is then i — -= Vaji^aj — C^jlb)- We emphasize again that Caj and 

C' • are assumed to be real in our calculation; i.e., CP noninvariance occurs only in the 



aj 



CKM matrix. If they are complex numbers, it will be completely another story |18 
The ordinary couplings in the Minimal Supersymmetric Standard Model ( MSSM ) are 
recovered by setting 

Caj = Ua cot P + dj tan j3, C'^j = —Ua cot P + dj tan /?, (3) 

where tan/5 is a parameter measuring the ratio of the vacuum expectation values of the 
two Higgs doublet fields. From now on, we always use the names of quarks to denote their 
masses. The vertices involving the would-be Goldstone bosons, UadjG^ and djUaG~, also 
arise as a special case: 

(^aj Ua Ctj, ^ aj ^ct Ojj. y±j 

To simplify the computation of diagrams involving W^ exchange, we shall use the back- 
ground field gauge |]T9| |^ [^ ( or the nonlinear R^ gauge [^ ) with ^ = 1. There will 



be no mixed W G^A vertex [ A is the external electromagnetic field ], and the Feynman 
rule for the W^W~A^ vertex is 

- ie[{kQ - k+)„g^p + {k+ - k^)^gp^ + (A;_ - h)pga^l], (5) 

where fco, ^+, ^- are incoming momenta for the fields A^, W'^ , ^a ■ Finally, to avoid 
any ambiguity, we define the effective EDM interaction as 

CeS = ^#75^^-^^^'. (6) 

where d is the EDM of the fermion ijj and it is real by Hermiticity. The Feynman rule for 
the effective vertex ipipA^ is 

-d[l^^4]l^^ (7) 

where q is the outgoing momentum of the photon. 
A. Renormalization of one loop elements 

We shall be mainly concerned with the EDM of the up-type quark. The result for the 
down-type quark will be obtained by simple substitutions at the end. To induce an EDM 
for the quark Me, the imaginary part of the CKM matrix must be involved so that the 
flavours j and k in Fig. 1 are different. Therefore, we need to renormalize the off-diagonal 
self-energy —iT}^{tj and the vertex with the photon ieV^M^, ^ + '?; <?)• Denote the bare one 
loop contribution by a hat and the counter-term one by a tilde, so that the renormalized 
quantities are 

Y^kj ^ fkj ^ fkj^ ^kj ^ fkj ^ fkj_ ^g^ 

The general structure of S^-' may be parametrized as 



t^^{i) = (\A{e) + B{e)^^\ + u^[c{e) + D{e)^^\. 



(9) 



The oflF-diagonal self-energy is renormalized by requiring that there be no mixing when 
either of d^ and dh is on-shell: 



d, 



H)Y!^\1)\^^^^=^, Y}^\l)dMi= 



0. 



(10) 



The counter-term is then determined to be 
1 



|]fci( 



d^ 



4) [d^A{d'^^ujO{d'M - [dkA{di)+u^C{di) 



dj^ 



1 



:ii^ 



dk + di 



dkhd-djB{d')+Uc,D{d')) + [dkB{di)+Uc,D{di))^,{^-d 



The renormahzation of the vertex T^^ is not independent but related to that of the self- 
energy by the Ward identity, 



q^teT'jii, i + q;q)= teQd {^"'{i) - ^"'{i + q) 



(12) 



It may be explicitly checked that the bare one loop quantities satisfy the above identity 
so that the latter must also be separately satisfied by the counter-term quantities. In this 
way, we find 



flK^,i + q-,q)=Qdl, 



+ 



dk + dj 



dk - dj 



(dkA{dl) - d,A{d])) + u^ {C{dl) - C{dl 



dkB{di) + djB{d')) + u^ [D{di) - D{d 



(13) 



75 



Since we shall present separate contributions from exchanges of H"^ and W^ in the sub- 
sequent subsections, we give below the functions A, B, C, D arising from exchanges of 
H"^^ W^ and G"^. We work in n = 4 — 2e dimensions to regularize ultraviolet divergences. 
For the H^ exchange. 



A{e) 
B{e) 

D{f) 



i4nr'Gp/V2V:,V^,{C^kC'^^ + C'^,C^j)h{e,ul,ml 
iAn)-^GF/V2V:,V^,{G^kCaj - C'^,C'^^)fo{f,ulml 
iAn)-^GF/V2V:,V^,iG.kC',^ - GUCaj)Mf,ulml 



(14) 



For the W"^ exchange in the ^ = 1 gauge, 

Aif) = -B{f) = +{An)-'g'v:,V^,\{2-n)f,{f,ulml.), 



(15) 



c{e) = D{e) = 0. 

The G^ exchange in ,^ = 1 gauge is a special case of the H"^ exchange, i.e., m'jj 

and with couplings C and G' substituted by values in Eqn.(4). The functions /i and /o 

are given in the Appendix. 



m 



w-> 



5 



B. Double H^ exchanges 

The momentum arrangement for external quarks and photon is shown in Fig.l. To 
pick out the EDM, we expand to the hnear order term in the photon momentum q. One 
should be careful in dropping terms that are superficially of zero order in q, since some of 
them are actually of linear order when the equation of motion is applied, and thus may 
contribute to the EDM. Notice that the final 75 in the effective EDM vertex can only 
come from Yukawa vertices since there would be no P violation if no 75 were involved 
in these vertices. To simplify the expression, the equation of motion for external quarks 
is freely used and only those terms that can finally contribute are kept. After a tedious 
computation, the Feynman diagrams in Fig. 1 sum up to the following structure with a 
common coefficient, {4:7Ty^eGp/2VekV*,.VajV*,q'''y5: 



+Ci{ea- kj) [Fi,^,(A;) - Fi,^,(j)] + C2{ea; kj) [F2,^.(A;) - F2,^.(j)] 

+C3{ea; kj)F3^^^{k) - C3{ea;jk)F3^^^{j) 

+C4(ea; kj) [F4,^^(A;) - ^4,^^(7)] + Cslea; kj) [F^^f^uik) 



5,fiu 



m , 



(16) 



where Ci are combinations of Yukawa couplings and quark masses dk and rf,-: 



Ci{ea;kj) 
C2{ea; kj) 

Csiea; kj) 

Csiea-Jk) 
C4,{ea; kj) 



[CekCej — C^j,C^,){C^^Caj — CakC^j) 



aj, 



- [yak'^aj — '^ak'^ajjyy'ek'^ej — '^ek'^ej 

Kyek^ei — ^ ek^ ei)y^ ak^ ctj " ^ak^^.^ 
+ [CakCaj — C^f^C^j){C^i^Cej 



did. 



jUk 



k i 



aj, 



Gp.k^^ 



-(CekCej — CekCej){Ce,kCaj — CakC^j) , , 



CpkC„ 



Y/ 

a; 

"'''dk 



4 



1 






[CakCaj CakCaj) [CekCej 

Csiea; kj)\j^k, 

— [CekCej — CekCej) [CakCaj + CakCaj) 



+ V^ek^ej — CekCej) [CakCaj + CakCaj)-!- 



k I "-J 
1 

1 



dk + dj 
1 



(17) 



dk, 



C^[ea]kj) = C4[ea;kj)\, 



dj_ 



The functions -Fj^^j, are complicated loop momentum integrals which are too lengthy to 
be displayed here. For brevity, we only indicate their dependence on the internal quark 
mass dk or dj although they depend as well on the external momentum p, the internal 
quark mass Ua and the charged Higgs mass itlh- 

The above results are obtained without using any approximations. The discussions 
in the next section will be based on these general results. We notice from Eqn.(17) that 
the coupling combinations Ci, C2, C4 and C5 are symmetric with respect to k and j 
such that the structure in Eqn.(16) is antisymmetric with respect to k and j. Therefore, 
the mirror-reflected diagrams corresponding to interchange of k and j in Fig.l are simply 



6 



related by V.kV^^V^.V^^ -> -KjV;* V^^V;* = - (KfcV^^fcV;,^*)*, so that in their sum the 
^^{^ekyak^oijyej) term is cancelled while the Im(V^fcV^afcV'ajV^*) term is doubled. This is 
essential, as emphasized for the SM case in Ref. |12[, to guarantee that the EDM is a 
real number as required by the Hermiticity of the effective action. The results for the 
down-type quark de are obtained by the following substitutions: 






Ve^/3*^ai^°e, 



r" 









Ur 



ak ^-' f3i^ "-" aj 

di, dj -^ Ua, dk -^ Uf3, Ue ^ 4- 



ej 



Before concluding this subsection, we would like to get some idea of how the quark 
EDM looks like. For this purpose, let us specialize to the case of the u and d quarks in 
the MSSM. We may use then the small external mass approximation ( SEMA ). In this 
approximation, only terms linear in the external mass are kept while higher order terms 
are safely ignored. [ At least one factor of external mass is involved due to the chirality 
flip feature of the EDM operator. ] The formula simplifies considerably. Fig. 1 along 
with its mirror reflection gives for the u quark. 



d{u) 
F{k) 



+eGllm{V^,V:,V^,V:MA7r)-Huul[F{k) - F(j)], 



+Qui 



+Qs 



d' 



dH 
(2^ 



^^(h{k)coep + h{k) 



- Ji_i cot^/3 



dl 



DhDI 



[h{k) cot^ (3 + Io{k)) + 



DHDk2 



(19) 



where Dh = i —mj^. For later discussion we should mention that the whole contribution 
from Fig. 1(a) is given by the term involving D]j^ upon setting Qu to 1. The /— and 
J— functions arise from the inner loop integration and their explicit forms are given in the 
Appendix. The relevant feature at the moment is that the J— functions depend only on 
£^ and the inner loop masses rrijj and u^ while the /—functions depend also on dl or d'j 
as indicated above. For the d quark. 



d{d) 
F{(3) 



+eGj,lm{V^,Vp,V:,V^,)iAn)-Hddi[FiP) - F{a)] 



+Qdi 



+Qui 



dH 

(2vr)4 

dH 


r 2 
Up 


{2^Y 


[DhD} 



(/i(/?)tan2/? + /o(/5) 



/i(/?)tan2/? + /o(/5)) + 



--Ji,i tan^/? 



DhDp2 

f /I 
I -J2 1 tan^ 3 + Jii 



(20) 



where D 



'P 



u 



73' 



Ji^i and J2,i depend on (?, d'f and m\ while /o(/?) and /i(/?) depend 
on m| as well. Suppose that m^, d^ = di and dj are respectively the top, bottom and down 
quarks, we find that, up to logarithms. 



d{u) ~ eGl6{4n) 



>mt 



nium^' 



m 



d{d) ~ eG'jpS{4:n) '^rridml 



m; 



m 



H 



^■(land cot2/5), 
( 1 and tan^ P 



(21) 



where S = lTsi{VubVflVtdV*^) . For numerical estimate we take the following input param- 



eters, Gf ~ 10-^ GeV-^ |(5| ~ 5 ■ 10- 
4.5 GeV, ruH ~ 200 GeV and tan/? 



m„ 



rrid 



5 MeV, mt ~ 170 GeV, m^ 



30, then. 



\d{u)\ -10^ 



-34 



e cm, 



\d{d)\ -10 



-31 



e cm. 



(22) 



Fo comparison, we quote the three loop result in the SM |jT5[, \d{u)\ ~ 0.35 ■ 10^'^^ e cm, 
\d{d)\ ~ 0.15-10"^^ e cm. Therefore, if not for the cancellation mechanism to be discussed 
in the next section, the result for the light quark EDM in the two Higgs doublet model 
would be very different from that in the SM. 
C. Outer loop W^ plus inner loop H^ exchanges 

Since we work in the background field gauge, we present the separate results from 
W^ and G^ exchanges. For the W^ exchange. Fig. 1(a) does not contribute to the 
EDM and Fig. 1(b) is completely cancelled by corresponding terms in diagrams (c)-(e). 
The counter-term diagrams do not contribute either. Diagrams (c)-(e) with their mirror 
reflection then give for the external Ue quark. 



+eg^V2GFlraiV,kV:,V^,V:^) (4 



vr) ^Mr 



1 . . fj^ 

2^ ^^J {2nyDw \Dk D, 






dk + dj 



^ak aj 



a] 



■"ak^OLj 



dk — dn 



TT ~ Tr) (Q^-Jo,i - QdJi,i) {^f,7u - 47m), 



(23) 



where Dw = {(- — v)"^ ~ i^w- '^^^ result for the external d^, quark is obtained by substi- 
tutions. In the SEMA, we have 



d{ue 



Veg^GF/V2{^'nrHm{Ve,V:,V^,V:^)\ 



ak aj 







Ub f\ (Jjp 



ak^aj 



D, 



dk — dj 

{QuJo,i — QdJl,l] 



(24) 



where now D^ = i'^ — m^^. The G^ contribution is a special case of the double H^ 
exchange; i.e., we only need to replace the couplings Cek, Cej, C'^^, C'^^ by their values in 
Eqn.(4) and Dfj by D]y. We display here the coupling combinations which are relevant 
to discussions in the next section: 



Ci{ea;kj) = 2uedkdj 



r^ r^ r^i r^i r^ r^i r^i r^ 

^ak^aj ^ak^aj ^ak^aj ^ak^oij 



dk + di 



dk - d^ 



C2{ea;kj) = —2u, 
C2,{ea]kj) = 2uedk 






■'ak^aj ^ ak^ 



aj 






dfc + dj KA,^ u,j 

■[{CakCaj — Cak^aj) ~ [CakC^j — C^j.Caj)\ 

Ciiea; kj) = -2ue{Cak - C'^k){Caj - C'^j), 

ul — dhdi 



(25) 



C^{ea]kj) 



Kyak^aj — ^ak^aj) . , . T V^ak^OLJ — ^ak^, 



, -Ug + dkdj 



"^ dk- d 



'3 J 



D. Outer loop H^ plus inner loop W^ exchanges 

The contribution from W^ exchange in Fig. 1 and its mirror reflection is 



+eg^V2GFlm{yekV:kVa,V:^){^n) 



^ek^ej — ^ek^ej ^ek^ej '~ ^ek'-'ej 



dk + di 



dk — dj 



(26) 



where Hfj,i, is another chain of loop momentum integrals. In the SEMA, we have for the 
Ue quark, 



d{ue) = +eg^V2GF{4iT)-^lm{VekV:,V^,V, 
1 



ep 



ek^ej — ^ek^ej ^ek^ej '~ ^ek'-^ej 



H{k) = +Quij 



dk + dj 
d'^i mird 



2 M 



dk — dj 



[H{k) - H{j)] 
dl 



" ' h{k) + Qs / 7K^,-fT^h{k) 



^" J {2nYD%Dk^''-^' ' ""V {27,YDhD^ 

1 . /• d'^i r 

[QdJ2,l — QuJl,l 



(27) 



+ 2'y {2nYDHDk 
f d* 

-{Qu - Qd)i j 



"•^2,0, 



{2nYDHD, 

where D^ = i"^ — m^j. The G^ contribution arises as a special case of the double H^ ex- 
change: m\ in the /— and J— functions is replaced by rn^ and the coupling combinations 
are given by the following ones, 



Ci{ea;kj) = —2uadkdj 






6*2(60;; A;j) = — 2-Uq, 
C3{ea;kj) = 2m«4 



dk + dj 



dk — fli 






ej 



ek'~-'ej 



dk + dj 



Ciiea; kj) 



iCkCe^-C'^kCe,^^^^^^ 



dk — dj 



I r^i \^k dj ^^ ^1 _ n' n \ "fc + dj 



[CekCej — C^j,C^j) 



ul - dkdj 
dk + dj 
C,{ea-kj) = -2u^{C,k-C'^k){Ce,-C'^j). 

E. Double W^ exchanges 



(2^ 



, (f^t /-i n ni \ ^" dkdj 



dk — d 



J J 



For completeness, we present in this subsection the result from double W^ exchanges, 
i.e., the SM result. Since we work in the background field gauge, we may separate four 
kinds of contributions: double W^ exchanges, double G^ exchanges, outer loop W^ plus 
inner loop G^ exchanges, and outer loop G^ plus inner loop W^ exchanges. For double 
W^ exchanges. Fig. 1(a) does not contribute and the contributions from (b)-(e) are 
completely cancelled. The counter-term diagrams do not contribute either. The case of 
double G^ exchanges is recovered from subsection B by mfj — > m^ and by evaluating 
couplings in terms of Eqn.(4); i.e., the coupling combinations become 

Ci(ea; kj) = 0, C2(ea; kj) = +8uaUe, C^iea; kj) = -SuaUedl, .^^. 

Ciiea; kj) = -Sw^Wg, C^iea; kj) = Suaul- 

The contribution from outer loop W^ plus inner loop G^ exchanges is obtained from 
Eqn.(23) by mfj -^ m^ and replacing the coupling combination in the square parentheses 
by — 4-Uq,. Similarly, for outer loop G^ plus inner loop W^ exchanges, we replace the 
coupling combination of Eqn.(26) by — 4-Ue- 

III. ANALYSIS OF CANCELLATION MECHANISM 

The SM result for the quark EDM was presented in the subsection 2E. For fixed 
internal quark fiavours a, j, k and external quark u^ and upon summing the pair of 
reflection-related diagrams, it has the following structure, 

dK) = Im(K.V;\V;,V;*) [H{k) - H{j)] , (30) 

where the two terms depend exclusively on flavours dk and dj respectively. [ Of course 
they also depend on the masses of quarks Wg, Ua and the exchanged bosons. ] When 
we evaluate the contribution involving exchange of charged Higgs bosons in the MSSM 
by using the couplings in Eqn.(3), we flnd that the above structure is also preserved. 
Summing over the three down-type flavours i, j, k while flxing the up-type flavour a, we 
arrive at 

rf(we) = +lm{V,kV:kV.,V:^) [H{k) - H{j)] + Im(K.K.* K..K1) \H{i) - H{k)] 

+Im(K,Kr,K..K1) \H{j) - Hii)] (31) 

= 0, 

where the second equality is due to unitarity of the CKM matrix; e.g., the H{k) term is 

"lm(V;fc\/4\4,K*) - Im(V;.\/^\4fc\/4)] H{k) 

:m [V^kV^ki^^e - V^kVX)] H{k) (32) 

= 0. 

The above cancellation occurs actually for any number of generations. However, it should 
be emphasized that the antisymmetric structure itself in Eqn.(16) and others does not 
guarantee the above cancellation. The crucial point is that the dependence on quark 

10 



flavours dk and dj is completely separate. We believe that this point is also responsible 
for the strong cancellation witnessed in the three loop QCD corrections in the SM. We also 
notice in passing that this cancellation is weaker than in the SM case where it occurs even 
before summation over flavours if one works in the unitarity gauge: Fig. 1(a) vanishes 
automatically and others cancel among themselves due to simple equalities [Q. In the 
present case however, Fig. 1(a) always contributes, e. g. as indicated in the subsection 2B, 
and there are no similar equalities which would demand the cancellation before summation 
over flavours is taken. Below we examine how this separate structure could be possible 
for general couplings Cai and C'^^. In other words, we want to determine what kinds of 
couplings are allowed for the separate structure to occur. 

Let us begin with the case of double H^ exchanges. It is natural to require from 
Eqn.(16) that Ci{ea;kj), C2{ea;kj), C^{ea-,kj) and C^{ea-,kj) be independent of d^ 
and dj and that C^^ea-jkj) can only depend on dk and C^^ea-jjk) only on dj. It is 
reasonable to assume that Cak and C'^f^ are universal as functions of quark masses d^ and 
Ua, and that these masses do not appear as denominators in Cak and C^^.. Then, we must 
have Ci{ea; kj) = 0, so that 

[Uek'^ej — ^ek^ejjyy'ak^aj — ^ak^aj) " Xyak^aj — ^ ak^ aj ) V^ ek^ ej ~ ^ek^ej)-i 
C2\ea; kj) = 2[CekCej — C^kC^j){CakCaj — CakC^j) ,2 _ j2' (oo\ 

Caiea; kj) = -4C2(ea; kj), 
Caiea-Jk) = Csiea; kj)\j^k- 

The independence in C2 of dk and dj along with assumptions about Caj and C' • implies, 

Cak(-^aj ~ C^k^aj ~ Va\dk + dj), ^ak^aj ~ ^ak^^j = Oa[dk — dj), (,'j4j 

where r]a and 6a are independent of dk and dj and may depend on «„. Furthermore, from 
the above equations, we have 

{Cak + C'ak){Caj - C'aj) = r]a{dk + dj) + 6a{dk - d^). (35) 

The factorized dependence on dk and dj on the left-hand side means that we may have 
two choices, 



Then, 



Case (I) : rja = +Sa, 

Cak + C'^k ^ dk, Cak - C'^k independent of 4; 
Case (II) : rja = —6a, 

Cak - C'ak oc dk, Cak + C'ak independent of dk 



C,{ea-kj) = 6e{CakTC'ak){Ca,TC'aj): 

C,{ea-kj) = 6a{CekTC',k)iCejTC',^), 



(36) 



(37) 



where the upper sign corresponds to the Case (I) and the lower sign to the Case (II). 
Therefore C4{ea; kj) and 6*5(60;; kj) are both independent of dk and dj. As one may 
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have realized, MSSM falls into the Case (I). Generally, the Case (I) corresponds to a left- 
handed theory in the sense that the Higgs bosons are doublets under SU{2)l, while the 
Case (II) corresponds to a right-handed theory under SU{2)r. As far as the contribution 
from double H^ exchanges is concerned, the two cases are equivalent up to a sign. This 
two-fold ambiguity is dismissed when contributions from mixed exchanges of H^ and W^ 
are considered, since W^ couples only to the left-handed current. Consider for example 
the case of the subsection 2C. The coupling combination in Eqn.(23) is {rja + Sa) while 
those in Eqn.(25) become 

Ci{ea;kj) = 2uedkdj{r]a - Sa), 

C-2{ea]kj) = -2ue{r]a + Sa), 

Csi^ea; kj) = 2uedk[r]a{dk + dj) + 5^(4 - d^)], (38) 

C^{ea-k3) = -2Me(afc-C;fc)(C«,-C^), 

C^{ea]kj) = 2[r]c,{ul-dkdj) + 6aiul + dkdj)]. 

The previous requirements on these combinations then single out the Case (I) as the 
only choice, otherwise the contribution in the subsection 2C cannot be cancelled upon 
summing over quark flavours and the quark EDM already arises at two loop order from 
mixed exchanges of H^ and W^. The subsection 2D does not give further constraints. 
The Case (I) may be parametrized as follows, 

^ak = Xadk + Ua, ^ak ~ ^adk ~ Vai ("j") 

where Xa and ya are independent of dk but may depend onua. A similar analysis may be 
repeated for the EDM of the down-type quark rfg, following the prescriptions in Eqn.(18). 
The cancellation of contributions from both double H"^ exchanges and mixed H"^ — W^ 
exchanges requires. 

Cat — C'ai '^ Ua, Cat + C'^i independent of Mq. (40) 

The constraints from Eqns.(39) and (40) then demand, 

Cai = xdi + yua, C'^i = xdi - yua, (41) 

where x and y are real constants that depend on the detail of the models and cannot be 
determined from this analysis. We see that in order to have a separate structure in the 
contribution to the quark EDM, the H"^ couplings with quarks cannot be arbitrary, but 
have to be in a form that is required by spontaneous symmetry breaking and the Yukawa 
couplings. 

IV. CONCLUSIONS 

We have presented an explicit two loop calculation of the quark EDM in the two Higgs 
doublet extension of the SM where CP noninvariance is assumed to be encoded in the 
nonzero phase of the CKM matrix. Naively we would expect a contribution which is much 
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larger than the SM result arising at three loop order. However, this large contribution is 
not realized in practice. Our detailed calculation shows that the contribution involving 
exchange of charged Higgs bosons vanishes completely due to a cancellation mechanism. 

We found that two factors are responsible for this complete cancellation. One is 
the unitarity of the CKM matrix, the other is the separate dependence in the relevant 
amplitude on the masses of internal quarks which are weak doublet partners of the external 
quark considered. We noticed that the antisymmetric structure itself of the amplitude 
in internal quark flavours is sufficient to guarantee the reality of the EDM, but is not 
to guarantee its complete cancellation. To analyse the cancellation mechanism in more 
detail, we examined the inverse problem of what kinds of couplings between iJ^ and 
quarks are allowed for the separate structure to occur. We found that the couplings must 
be in a form that is dictated by spontaneous symmetry breaking and the Yukawa coupling. 
This is consistent with our assumption that the origin of the CP noninvariance resides in 
the complex CKM matrix. 

/,From the vanishing result of the quark EDM at the lowest nontrivial order, it is safe 
to conclude that it should be difficult to detect the CP noninvariance through the EDM 
of the neutron or the leptons if the CP noninvariance is of the CKM origin. 
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APPENDIX A: SOME FUNCTIONS APPEARING IN ONE LOOP DIAGRAMS 

The functions fi {i = 0, 1) appearing in renormalization of one loop elements are 

h{i^, m^, M^; x) = xM'^ + (1 - x)m'^ - x{l - x)f, 

where m and M are respectively the masses of the quark and the boson in the loop. The 
on-shell subtraction then produces the functions Iq and Ji which depend on dl or dp as 
well as on £^, m'^,M'^. These functions generally involve the function h to the power of 
— e; but in the SEMA, it is sufficient to expand them to the zero order in e (i = 0, 1): 



The functions Jij{i = 0, 1, 2; j = 0, 1) arise from differentiation with respect to £^, 
or M^ of the functions fi'. 



2 

m 



J,, = J,,{f,m\M')= fdx-p^^^^ + 0{e). (A3) 
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Figure Captions 



Fig. 1 Diagrams that contribute to the EDM of the up-type quark Ue- The wavy lines 
represent the electromagnetic fields and the dashed lines the H^ or W^ fields. Diagrams 
for the down- type quark de are similar, with the replacements: a ^ i and j, k ^ a, (3. 
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